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Abstract 

The asymptotic shape of randomly growing radial clusters is studied. We pose the problem in 
terms of the dynamics of stochastic partial differential equations. We concentrate on the properties 
of the realizations of the stochastic growth process and in particular on the interface fluctuations. 
Our goal is unveiling under which conditions the developing radial cluster asymptotically weakly 
converges to the concentrically propagating spherically symmetric profile or either to a symmetry 
breaking shape. We demonstrate that the long range correlations of the surface fluctuations obey 
a self-affine scaling and that scale invariance is achieved by means of the introduction of three 
critical exponents. These are able to characterize the large scale dynamics and to describe those 
regimes dominated by system size evolution. The connection of these results with mathematical 
morphogenetic problems is also outlined. 

PACS numbers: 05.40.-a, 05.65.+b, 87.10.-c, 87.10.Mn 



I. INTRODUCTION 



The building blocks of mathematical morphogenesis were put several decades ago in the 
seminal works of Thompson pQ and Turing [2J. A particularly relevant problem in this 
context is the examination of the properties relating to the architecture of cell colonies, as 
already noted and investigated by Eden [3, H] . At some basic underlying level, one could 
say that Turing and Eden shared the goal of understanding how a macroscopic structure, in 
particular one breaking the initial homogeneity, could arise out of a multiplicity of simple 
interactions. In the introduction of [2] one can read Turing thoughts on how a structure 
composed by randomly and isotropically proliferating cells could or could not break the initial 
spherical symmetry. But then his approach moved to the use of reaction-diffusion equations 
and to the search of a deterministic mechanism, the currently well known Turing instability, 
able to give rise to pattern formation. On the other hand, Eden focused his interest on 
the evolution of a pure growth process. He concentrated on a probabilistic abstraction of 
a developing cell colony and studied stochastic symmetry breaking concomitant to growth. 
In particular, Eden studied the architecture of a lattice cell colony to which new cells were 
added following certain probabilistic rules. The objective was determining the asymptotic 
colony profile. If the growth rules are isotropic then the realizations of the growth process 
are either spherically symmetric or this spherical symmetry is broken by means of long range 
fluctuations. 

These sorts of discrete models are usually studied using continuum equations for ease 
of analytical treatment (5j E]. While most of these developments have been carried out 
for systems in which there is an external input of mass, the universality of the surface 
fluctuations is expected to remain unchanged when the origin of mass appearance is internal. 
In fact, the Eden model in a strip geometry has been satisfactorily analyzed within the 
classical theoretical framework [7j. And, as we already have mentioned, this model was 
introduced as a probabilistic abstraction of a proliferating cell colony. So we expect our 
results to hold independently of the internal or external origin of mass input into the system. 
This will be so at least in those cases in which the growing radial interface can be described as 
a Monge patch in spherical coordinates. Growing clusters as the ones generated by diffusion- 
limited aggregation processes yield more complex structures that cannot be described in such 
simple terms, and thus fall beyond the scope of the present approach. 
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Indeed, the original Eden problem can be greatly generalized by means of the use of 
stochastic partial differential equations. They allow a systematic study of the properties 
of the colony periphery, particularly of the interface fluctuations. At the same time, they 
allow us to get rid of the undesirable lattice anisotropy. In this work we will concentrate 
on the properties of the realizations of the stochastic growth process. Our goal is unveiling 
under which conditions the developing radial cluster asymptotically weakly converges to the 
concentrically propagating spherically symmetric profile. Let us emphasize that we are in- 
terested in determining the presence or absence of spherical symmetry in every realization of 
the stochastic growth process in the long time limit. As we are considering isotropic growth 
which is free from deterministic instabilities, averaging the cluster profile over many real- 
izations yields an immediate spherical symmetry. The properties of the realizations can be 
ascertained, as we will see, calculating suitable correlations adapted from classical elements 
of stochastic growth theory [7j, and keeping in mind that we should expect weak rather 
than strong convergence of the asymptotic profiles. Weak convergence is a mathematically 
totally precise concept which physically corresponds to the convergence resulting from the 
self-averaging of wildly oscillating quantities or structures [8] . 

Apart from the overall shape, the microscopic fluctuations of growing radial clusters have 
been studied as well. It turns out that these fluctuations sculpture a fractal surface which 
is statistically self-affine [7J |9j [TO]. This microscopic roughness, in those cases in which 
the pointwise width [TT] is well defined (which actually rules out all flat interfaces due to 
the microscopic properties of white noise) is described by the familiar scaling of planar 
systems pJ2]. This is not so, however, for long range fluctuations p2]. Long range radial 
correlations cannot be in general deduced from planar scalings. They still adopt a self-affine 
form for macroscopic length scales, but scale invariance is this time achieved for different 
values of the critical exponents. The following sections will discuss the numerical values and 
physical meaning of these exponents. We will see that they are able to characterize which 
regimes asymptotically weakly converge to the radially symmetric shape and which do not. 

II. RADIAL CORRELATIONS 

In planar situations of non-equilibrium growth the interface is parameterized by the height 
function h(x,t), which indicates the displacement of the interface with respect to some 
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hyperplane taken as origin, at some d— dimensional spatial point x and for some instant of 
time t [7]. We illustrate this geometry in Fig. 0, where we present the one-dimensional 
solution to the so-called random deposition equation 

d t h = F + t(x,t), (1) 

where £ is a zero-mean Gaussian noise whose correlation is (£(x, t)£(x', t)) = e5(x—x')5(t—t'). 
This numerical approximation is found by means of spatially discretizing the system. 
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FIG. 1: Solution to the random deposition equation (JIJ at several times. The values of the 
parameters are F = 4, e = 1, and the system has been discretized in 2 10 spatial points. We have 
used as initial condition h(x, 0) = 0. 



In these systems and for large spatiotemporal scales the height-height correlation function 
adopts the self-affine scaling form 

(h{x,t)h{x',t)) = t Wc f fu '~ r 

= \x- x'\ 2a °g {^j^ 1 J , (2) 

which is invariant to the transformation x — > bx, t — » b z t and h — > b a h, for a = (3z and 
a = a c , (3 = f3 c and z = z c , where the subscript c denotes the classical value of an exponent 



for a given model, a is known as roughness exponent, (3 as growth exponent and z as 
dynamic exponent, / and g are scaling functions. This scaling is valid for a large class of 
linear and nonlinear models. If we consider the linear Langevin equation 

dth = -D\V\ c h + £(x,t), (3) 

where £, as before, is a zero-mean Gaussian noise whose correlation is (£,(x,t)£(x',t)) = 
eS(x — x')5(t — t') and | V|^ is a pseudo-differential operator to be understood in the Fourier 
transform sense, the exponents read z = (, a = (£ — d)/2 and /3 = 1/2 — d/(2(). One 
would like to know if the radial counterpart of Eq. ^ displays an analogous behavior to 
that dictated by scaling (|2]). In the radial extension of the theory one has a new degree of 
freedom which was trivial in the planar case: the rate of growth. If we assume that the 
average radius (measured as distance from the origin) of the macroscopic radial form grows 
as a power law (r) = Ft 1 for t > t and a growth index 7 > 0, then the equation for the 
radial surface fluctuation p(9, t) is 

„ D c 1 r]{6,t) 

according to the reparametrization invariance principle [HI [T5|, and where 9 denotes the 
set of angles parameterizing the d— dimensional radial interface and J (6) is the Jacobian 
determinant of the change of variables from (x,h) to (0,r) evaluated at r = 1 |16j . The 
noise is simply 77 = ^/v^e. In order to visually compare planar and radial processes we have 
plotted the field r(6,t) = (r) + y/ep(6,t) for D = in Fig. Q. Neglecting the diffusion 
constant makes this stochastic process the radial counterpart of the random deposition 
process plotted in Fig. (fil). We will therefore refer to it as the radial random deposition 
process. 

As we have mentioned in the introduction, one of the motivations of this work is studying 
the large scale, long time properties of a randomly proliferating set of cells or entities in 
general. Such a process can be modelled in some detail using a master equation, which can 
be mapped into a stochastic partial differential equation provided with a multiplicative noise 
term [T71 HE]. On the other hand, we focus in this work on stochastic partial differential 
equations provided with additive noise. While this could seem as an apparent paradox, the 
underlying reason is that equations ^ and Q are effective models. So they are intended to 
describe just the long range properties of the growing interface, and in this limit fluctuations 



FIG. 2: Radial random deposition growth process plotted at different times. The values of the 
parameters are F = 2, 7 = 1, e = (2007r) _1 , and the system has been discretized in 2 11 spatial 
points. We have used as initial condition r(9, 1) = 2 (note that we have used as initial time to = 1). 



can be effectively considered additive. A nice example of this, a case in which a front driven 
by multiplicative noise is found to belong to a universality class described by an additive 
noise stochastic equation, is described in ] ] 9] . 

The analysis of the exact solutions to Eq. Q for 7 > l/£ and long times shows the 
diffusion term is irrelevant in the large scale \0 — 6'\ t 1 ^~' y and the correlation becomes [To! 



C(9,9',t) ~ < 



where C(9,9',t) 



t 1 ~ ld 5{9 — 9') JJ{9) if 7d<l, 
\n(t)S(9 - 9')/J(9) if 7 d = 1, 
5{9-9')/J{9) if 7 d>l, 

(p(9,t)p(9' ,£)). In terms of arc-length differences i 



(5) 



P{9 - 9') one 
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if 7<i < 1, 
if jd = 1, 
if jd > 1. 

',£—>■ and p — >■ 6 Q p for z 



(6) 



finds [16] 

t5{£-£')/J{t^£) 
C(£,£',t)~ I t\n(t)5(e-£')/J(t-<£) 

P d 5{£-£')/J{t-U) 
This correlation is invariant to the transformation £ — >■ 
I/7 =: z r and a = 1/(27) — c?/2 =: a r (a = =: a r ) when 7c/ < 1 (jd > 1). We define 
/3 r := a/z = a r /z r = (1 — r yd)/2 (f3 r = 0) when 7c? < 1 (7c? > 1). First, the functional form 
of this correlation is not reducible to scaling form ([2]), and second, the exponents are totally 
different (this makes precise the analysis carried out in [13] ) . In summary, long range radial 
surface fluctuations are not describable in terms of scaling (|2]). 



III. DILUTION AND UNIVERSALITY CLASS BIFURCATION 

It is clarifying to consider a simplified problem: an abstraction in which Eq. ^ is placed 
in a growing domain [20] • If this is realized by means of the transformation x — > [t/t^x 
one finds 

dth = -D (I) ^ \V\ c h + 7FF- 1 + (I) d1 ' 2 £, (7) 

where the mass source is explicitly included and we assume fast growth 7 > l/( and the 
rough interface inequality ( > d. In this case it can be shown that the correlation obeys 
scaling u2« for short length scales \x — x'\ <C ^( 1_ C7)/C anc i f or long length scales \x — x'\ ^> 



f(l-C7)/C it is m 

' tS(y - y') if !d < 1, 

C{y,y\t) ~ \ t\n(t)5(y - y') if 7 d=l, (8) 
P d 5(y - y') if jd > 1, 

where y = (t/t^x is the counterpart of the arc-length coordinate. This correlation is scale 
invariant for the exponents {a r , z r } but not for the exponents {a c ,z c }. 

This way of realizing domain growth agrees with reparametrization invariance [16] , but is 
not unique. It makes space and mass grow simultaneously [16], while it is possible to realize 
domain growth at constant mass by introducing a dilution term [2HI EI] 

dth = -D (I) ^ |V| c /i -*lh + !Ft^ + (|) ^ e- (9) 
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In this case the correlation is 

C(y,y',t)~t5(y-y'), (10) 

which is scale invariant for both {a c , z c } and {a r = 1/(27) ~~ d/2,z r }. The first invariance 
can be used to argue that this correlation is a direct consequence of scaling (|2]), but the 



scaling form (10) is definitely ambiguous in view of the second invariance. Both sort of 
dynamics, with and without dilution, are physically realizable; while dilution dynamics is 
apparently better suited to describe some natural processes, no-dilution dynamics seems 
to agree better with the numerical simulation of some discrete models as discussed in [16]. 
Together with the difference in the roughness exponent, very pronounced differences in the 
auto-correlation and persistence exponents separate the interface behavior in presence and 
absence of dilution [20J. This reveals that the non-uniqueness of the extension of a Langevin 



equation to a growing domain implies in turn a universality class bifurcation. 

IV. GEOMETRIC CONSTRAINTS 

One can apply what has been learned in the previous section about growing domains to 
the radial situation. Including dilution in the radial equation for interface motion yields a 
different dynamics from the one established by reparametrization invariance. The correlation 
in this case becomes [H] 

which is scale invariant for the exponents {a r ,z r } but not for {a c ,z c }. The difference 
among a r and a r originates in the memory effects that affect the surface described by 
reparametrization invariance and disappear by virtue of dilution [161 EOJ [22] . The difference 
between the sets {a r , z r } and {a r , z r }, and the classical set {a c , z c } is of a different nature. 
Note that the Langevin dynamics on a growing domain and influenced by dilution gave rise 



to the random deposition correlation (10): 



(h(y,t)h(y',t)) ~t5(y- y'). (12) 
This correlation is invariant to the transformation 

y -»• by, t -> b z t, and h -»■ b a h, (13) 
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for 

z = k and a—(K — d)/2, kgR; (14) 

let us remark that scale invariance is achieved for any real number k, thus (13) and (14) 
actually constitute a one-parameter family of scale transformations. This is, there is one 
degree of freedom in the critical exponents, or in other words, the scaling form is insensitive 



to the dynamic exponent. Contrarily, in a radial geometry the correlation, Eq. (11), is 
scale invariant only if k = I/7: this is the constraint imposed by the geometrical factor 
J. The Jacobian determinant J = sin(6'i) d ~ 1 sin(^2) a! ~ 2 • • • sin(#<i_i) for a d— dimensional 
interface, where #1, • • • , Qd-\ is the set of polar angles, while it is of course independent 
of the azimuth angle So this constraint appears for d > 2, but it is not present for 
d = 1. The reason is that the one-dimensional radial interface is topologically identical to 
the one-dimensional, topologically toroidal, interface which results from assuming periodic 
boundary conditions on the growing domain setting. In this case, the angle-like coordinate 
x, to which the azimuth angle is akin, may be either dimensional or dimensionless (this 
could be expressed by means of the dichotomy of Lagrangian/Eulerian coordinates similarly 
to the developments in [SU] )• On the other hand, the polar angles, ranging only it radians, 
cannot be associated to periodic coordinates which could have dimensions. This fixes the 
scaling z = I/7. 



One can use the ambiguity of Eq. ( 12 ) to claim that z = ( is a good dynamic exponent 
in this case, but that would be accepting that topology can modify the scaling. Also, 
the value of the auto-correlation exponent in this case suggests that z = I/7 is the good 
dynamic exponent [20J. In any case, it is now clear that there is no basis for claiming 
that correlation ( 10 )-( 12 ) is a consequence of scaling Q in the present framework. This is 



because, although both forms are invariant to the transformation (13) for z = ( (for linear 



equations like (p|), self-similar form ( 10 )-( 12 ) preserves scale invariance independently of 



the value of z (as specified in Eq. (|14[)). So this second form is much more general and the 
value z = C plays no special role in it (let us note again that the value 2 = 1/7 plays a 
special role in the sense that it appears in the auto-correlation function [2D])- Despite this 



ambiguity in the growing domain situation, it is clear that the radial correlation (11 ) is just 



invariant for {a r ,z r }, and the classical exponents {a c ,z c } play no role in this case. This 



can be seen by means of a direct application of transformation (13) to this correlation and 
substituting a and z for the two different alternatives. 



V. DIMENSIONAL ANALYSIS AND NONLINEAR DYNAMICS 



Our results so far allow us to develop a consistent dimensional analysis of radial equations. 
We start with Eq. ^ performing the scale transformation p — > b a p and t — > b z t. The angular 
variables are dimensionless and this, as specified by the factor J, fixes the exponent z — I/7. 
The stochastic term becomes scale invariant for a = 1/(27) — d/2 and the diffusion irrelevant 
for 7 > 1/Q: this agrees with our previous analysis using the explicit solutions of the linear 
equations. We can now use this technique to get some insight on nonlinear equations. The 
Kardar-Parisi-Zhang (KPZ) equation [23J for the radial surface fluctuation can be written 
as 

d tP=^ d eP + ^(de P ) + Fd/2pd/2 J(ey/2 . (15) 

We note that scaling arguments were successfully employed in determining the critical be- 
havior of KPZ with a time dependent coefficient of the nonlinearity [23]. If we select the 
same values of a and z as in the linear case in order to make the stochastic term scale 
invariant we find that diffusion and nonlinearity are irrelevant for 7 > 1/2 and 7 > 3/(4 + <i) 
respectively. This agrees with the naif dimensional analysis of KPZ and thus we cannot 
expect it to be correct. In fact, the naif dynamical exponent z = (4 + d)/3 found by means 
of a dimensional analysis of the KPZ nonlinearity does not agree with the one measured 
in simulations and successfully predicted by the dynamic renormalization group analysis in 
one dimension [7]. It is plausible that the nonlinearity interacts with the noise to yield a 
correlation different from ^ (or from (11) if dilution is included). In order to guess what 



form the correlation would have we can try to get some illumination from scaling ([2]). Taking 
the short time limit in ^ we find 

C(x, x', t) ~ t 2 ^ +d/z <5(x - x'), (16) 

where the temporal prefactor expresses the time dependence of the variance of the interface 



center mass position [16] . Using this form we modify, for instance, Eq. (11) to the new 
correlation 



c'(£,£^t)~t 2 ^/-^ F ^, a?) 

which yields the exponents a = (2(3 C + d/z c )/(2j) — d/2, 2 = 1/7 an d A- = (2/3 c + d/z c )/2 — 
d^y/2. In this case we have the diffusion is irrelevant whenever 7 > 1/2 and the nonlinearity 
is irrelevant if 7 > l/z c , where we have employed the relation a c + z c = 2 in the derivation of 
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the last inequality. This last inequality could be considered as desirable but correlation (17) 
has a different drawback. As we have said, the exponent /3 r describes the time power law 
dependence of the standard deviation of the interface center of mass position. For the Eden 
model, which in planar format belongs to the KPZ universality class [7J, d — 7 = 1 and so 
P r = 1/6. This exponent however does not agree with the one measured in simulations [25] 
and this constitutes another disadvantage of the use of scaling ^ in radial systems. We 
note that considering the planar KPZ equation on a growing domain does not simplify much 
things with respect to considering the full radial case [25] . 

Together with the KPZ equation, the Villain-Lai-Das Sarma (VLDS) equation [271 EE] 
is one of most important nonlinear models for surface growth. If one derives the VLDS 
equation in the radial setting one finds irrelevance of the diffusion for 7 > 1/4 and of the 
nonlinearity for 7 > l/z c - In this case the well known hyperscaling relation for conserved 



growth models yields 2/3 c + d/ z c = 1, and so equality of (17) with ( 11 ). This suggests that the 
scaling properties derived herein for linear models are also valid for some nonlinear models 
as the VLDS equation. 

VI. CONCLUSIONS 

In this work we have examined the dynamic scaling of radial interfaces. We have shown 
that long range fluctuations, which are describable by two-point correlation functions, be- 
come scale invariant for values of the critical exponents which may totally differ from the 
corresponding ones in a planar geometry. In rapidly growing radial systems we have found 

Id 1 1 7^ , . 

a — , z = — , p r = , (18) 

2 7 2' 7' H 2 2 ' v ; 

all of which have a clear physical meaning. The j3 r exponent, as we have discussed, describes 
the interface center of mass fluctuations. The a exponent associated to the radial surface 
fluctuation indicates how much the growing cluster deviates from a d— dimensional growing 
hypersphere. For a < the growing cluster converges to a hypersphere concentrically 
growing at the deterministic rate t 7 , and in this sense the growth profile is flat. This sort 
of convergence is of course weak: the profile converges to the hypersphere on (any type 
of spatial) average. The microscopic roughness, which is given by the height-difference 
correlation, reduces to that of the planar profile [20] and precludes strong or pointwise 
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convergence. This is so at least in situations in which the height-difference correlation is 
well defined; otherwise a different scaling is possible even locally in space [TO]- This makes 
precise and generalizes our previous results [T3] in the following way: a fast growth rate will 
make the interface weakly converge to a hypersphere. We have further found the scaling 
relation a = (3 r z, which quantitatively relates the center of mass fluctuations, growth rate 
and rate of (weak) convergence to/divergence from the hyperspherical profile. These results, 
together with our previous ones (161 [20] , suggest the following generalization of the critical 
exponents for arbitrary values of 7 

«=^, V-=- = \-^, z = minkl|, (19) 

2 z 2 2z I 7 J 

which should be valid for linear models as well as some nonlinear equations as the VLDS one, 

which present a conserved drift and its critical behavior is accessible to dimensional analysis. 

For the KPZ equation, which nonlinearity introduces both mass and fluctuations into the 

system, we expect a different behavior. As we have seen, the amount of mass entering the 

interface (and of course the way it is introduced) may change the critical behavior. We have 

explicitly shown this by means of the linear models in absence and presence of dilution, 

which gave respectively rise to correlations ^ and (11). 



7<i 


h <7< a 


i <7< dh 


7 >max{l,^} 


Classical exponents 


Roughness 


Super-roughness 


Weak convergence to the hyperspherical profile 



TABLE I: Morphologies for different spatial dimensions and exponents. 

We have summarized the different morphologies that can appear in table |T] In the first 
column, for small 7, we have placed the regime in which the classical values of the exponents 
are recovered. The other three columns represent situations in which this is not the case. 
The "roughness" column describes the situation in which the radial fluctuations do not 
average each other out for long times, so in this regime the overall appearance of the cluster 
deviates from the hyperspherical profile. The "super-roughness" column states a particular 
case of this last regime, in which the amplitude of the radial fluctuations grows faster than 
the arc-length sustaining them. Finally, the last column describes the regime in which the 
weak convergence to the hyperspherical profile is reached asymptotically in time. 



Note that the meaning of the exponent f3 in (19) changes for 7 < l/z c and for 7 > 



l/z c . For 7 < l/z c this exponent describes the short time dependence of the spatially 
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averaged pointwise width with time, as in the classical situation. For 7 > l/z c , as we have 
mentioned, it describes the center of mass fluctuations, i. e., the less dominant source of 
fluctuations is chosen in each regime. For 7 < l/z c and a > (which implies the height 
difference correlation is well defined and the pointwise width is finite) there are always 
center of mass fluctuations. For 7 > l/z c the situation is different as it is possible to find 
growth regimes with no center of mass fluctuations and which pointwise width increases as 
a power law of time. This is precisely the weak converge to the hyperspherical profile: while 
microscopically the interface is rough, macroscopically it is a concentrically propagating 
hypersphere. The change of meaning of the exponent /3 coincides with a change of physics. 
The long range shape of the interface is no longer affected by the diffusion mechanism, 
which becomes irrelevant, but for the system size evolution, which sculptures the resulting 
profile [20]. And so, the violation of scaling ^ is associated with the advent of a new 
dominant physical mechanism in the interface macroscale: system size growth overtakes 
diffusion and determines the cluster macroscopic shape. 
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